Rough Path Analysis Via Fractional 

Calculus 

Yaozhong Hu* and David Nualart^ 
Department of Mathematics , University of Kansas 
405 Snow Hall, Lawrence, Kansas 66045-2142 



Abstract 

Using fractional calculus we define integrals of the form f{xt)dyt, 
where x and y are vector-valued Holder continuous functions of order 

/?£(-, — ) and / is a continuously differentiable function such that /' 

is A-H61dr continuous for some A > ^ — 2. Under some further smooth 
conditions on / the integral is a continuous functional of x, y, and the 
tensor product x ® y with respect to the Holder norms. We derive some 
estimates for these integrals and we solve differential equations driven by 
the function y. We discuss some applications to stochastic integrals and 
stochastic differential equations. 



1 Introduction 

The theory of rough path analysis has been developed from the seminal paper 

by Lyons 4 . The purpose of this theory is to analyze dynamical systems 

dxt = f{xt)dyt, where the control function y is not differentiable. If the rough 

control y has finite p-variation on bounded intervals, where p > 2, then the 

dynamical system is a continuous function, in the p-variation norm, of y and 

fc 

the associated multiplicative functionals y (E> ■ ■ ■ y, with k — 2, . . . ,[p\. In the 
case I < p < 2, the dynamical system can be formulated using Riemann-Stieltjes 
integrals and applying the results of Young In this case, Xt is a continuous 
function of y in the p- variation norm (see Lyons jS]). 

Suppose that / and g are Holder continuous functions on the interval [a,b], of 
order A and /i, respectively, with X+fj, > 1. Then, the Riemann-Stieltjes integral 
fdg can be expressed as a Lebesgue integral using fractional derivatives (see 
Zahle in] and Proposition 2.1 below). This fact has been exploited by Nualart 
and Ra§canu in ^ to analyze dynamical systems driven by a control function 
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y which is Holder continuous of order /3 > i. In this case further results are 
obtained in |2 along the line of the present paper. 

The purpose of this paper is to analyze dynamical systems dxt — f{xt)dyt, 
where the control function y is Holder continuous of order /3 G (ji 5)1 using the 
techniques of the classical fractional calculus, and following an approach inspired 
by the work of Nualart and Ra§canu ^ in the case (3 > \. In order to achieve 
this objective, we first provide in Section 3 an explicit formula for integrals of the 
form f{xt)dyt, where x and y are Holder continuous of order f3 G (|j 5)- This 
formula, given in Theorem 3.1, is based on the fractional integration by parts 
formula, and it involves the functions x, y, and the quadratic multiplicative 
functional x ® y . Notice that this explicit formula does not depend on any 
approximation scheme. As a consequence, we derive estimates in the Holder 
norm for the indefinite integral. 

Section 4 is devoted to establish the existence and uniqueness of a solution 
for the dynamical system dxt = f{xt)dyt. The main ingredient in the proof of 
these results is to transform this equation into a system of integral equations 
for X and x ® y that can be solved by a standard application of a fixed point 
argument. We show how the solution depends continuously on the Holder norm 
of y and y®y. We also prove some stability results for the differential equations 
which are interesting, new and may be difficult to obtain by other approaches. 
Remark that to derive our results we do not make use of the theory of rough 
paths, and we obtain explicit formulas that do not depend on any approximation 
argument. 

These results can be applied to implement a path-wise approach to define 
stochastic integrals and to solve stochastic differential equations driven by a 
multidimensional Brownian motion. As an application of the deterministic re- 
sults obtained for dynamical systems we derive a sharp rate of almost sure 
convergence of the Wong-Zakai approximation for multidimensional diffusion 
processes. We couldn't find this kind of estimates elsewhere. Similar results hold 
in the case of a fractional Brownian motion with Hurst parameter H S (5, 
The approximation of the solutions of stochastic differential equations driven by 
a fractional Brownian motion with Hurst parameter G (i, i) is more involved 
and it will be treated in a forthcoming paper. 



2 Fractional Integrals and Derivatives 

Let a, 6 G M with a < b. Denote by {a,b), p > 1, the space of Lebesgue 
measurable functions / : [a, 6] — > M for which ||/||^p(^ < 00, where 

[ ess sup {1/ : t e [a, 6]} , if p = oo. 
Let f ^ {a,b) and a > 0. The left-sided and right-sided fractional Riemann- 
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Liouville integrals of / of order a are defined for almost all x G (a, &) by 
and 



respectively, where (—1)^" = g-iu-a g^^^^^ T (a) = r"~^e~^dr is the Euler 
gamma function. Let I"^{LP) (resp. Ij^_{LP)) be the image of LP{a,b) by the 
operator (resp. If / e (L^) (resp. / G (L^)) and < a < 1 

then the Weyl derivatives are defined as 



^ ' r(l-a) I (i-a)" A (i-s)"+i ' ^ ^ 



and 

n" (-1)" ( fit) f'm^IlEl,] (00^ 



where a < t < b (the convergence of the integrals at the singularity s = t 
holds point- wise for almost all i G (a. fe) if p = 1 and moreover in L^-sense if 
1 < p < cxd). 

For any A G (0, 1), we denote by C^ia,b) the space of A-H61der continuous 
functions on the interval [a, b\. Recall from [7] that we have: 

• If a < i and q = -r-^ — then 

p ^ 1 — ap 



4\ (Lf ) = (if ) C (a, 5) 



• If a > i then 

p 



4\(L'') u/^(L^')cC""Ma,fe)- 

The following inversion formulas hold: 

la+iK+f) = /, yf^I^+iLP) (2.3) 
Ia-{DZ^f) = f, Vfei:_iLP) (2.4) 

and 

D:+{i:+f)^f, D:_{i:_f)^f, V/ Gil (a, 6). (2.5) 
On the other hand, for any /, g G L^ia, b) we have 

i:+fit)9it)dt = (-1)" fit)I^_git)dt , (2.6) 
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and for / G 1^+ {LP) and g £ (LP) we have 

D^+f{t)g{t)dt - (-1)-" f{t)Dt_g{t)dt. (2.7) 



Suppose that / e C^{a,b) and g £ C^{a,b) with A + /.t > 1. Then, from 
the classical paper by Young [5], the Riemann-Stieltjes integral fdg exists. 
The following proposition can be regarded as a fractional integration by parts 
formula, and provides an explicit expression for the integral fdg in terms of 
fractional derivatives (see 



Proposition 2.1 Suppose that f e C^{a,b) and g G CP^{a,b) with X + fi > 1. 
Let X > a and fi > 1 ~ 
it can be expressed as 



Let X > a and > 1 — a. Then the Riemann Stieltjes integral J^^ fdg exists and 



fdg = i-ir I D2,f (t) Dlz"gb- (t) dt, (2.8) 



where gf,^ (t) = g (t) - g (b). 



We will make use of the following two- variable fractional integration by parts 
formula, whose proof is given in the Appendix. 

Lemma 2.2 Let 'fi{S,,ri) and ip{£,,ri) be two functions of class defined on 
o- < ^ < r] < b. Suppose ip{^,ri) vanishes on the diagonal. The following 
fractional integration by parts formula holds for any < a < 1 . 

I' d^ I' V)didg = - £ DlfD'^^J^a+,b- (C, v)dL 

where D^^ denotes the fractional derivative on variable D'^'^ denotes the 
fractional derivative on the variable rj, and the operator F" is defined by 

V^rPiLv) = A':"'^^4+"'VXC,^)- (2.10) 

3 Integration of Rough Functions 

Fix i < < i. Suppose that x : [0, T] ^ M™ and y : [0, T] R"^ are /3-Holder 
continuous functions. Following [l] we assume that x (g)y is well-defined and it 
is a continuous function defined on A := {(s, t) : < s < t < T} with values on 
M™ (g) verifying the following properties: 

i) For all s < M < i we have (multiplicative property) 

® y)s,u + (a; ® y)u,t + - Xs) (S> (yt - y«) = {x ® y), ^ . (3.1) 
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ii) For all (s, t) e A 

{x<E>y),^t <fc|t-s|2^. (3.2) 



That is, {x,y,x (E) y) constitutes a multiplicative functional in the sense of 
the rough paths analysis theory. We will say that {x,y,x (8> y) is a (3-Hdlder 
continuous multiplicative functional on M™ OM''. 

If X and y are smooth functions, then 

{x®y)%= [ dx\df^ (3.3) 

clearly defines a /?-H61der continuous multiplicative functional. 

Let / : R™^ R'^ be a continuously differentiable function such that /' is 
A- Holder continuous, where A > ^ — 2. Our aim is to define the integral 

/ f{xr)dyr^Y. / hMdyl (3.4) 

using fractional calculus. 

Fix a number a such that 1 — /3 < a < 2/3 and a < This is possible 

because 3/3 > 1 and > 1 - /3. 

Notice first that the fractional integration by parts formula H2.8|l cannot be 
used to define the integral (|3.4|) because the fractional derivative -0"+/ (x) is not 
well-defined under our hypotheses. For this reason we introduce the following 
compensated fractional derivative: 



1 (1 — a) \ (r — a)" 

L J ■ ^ 

This derivative is well-defined under our hypotheses because 

1/ {Xr) - f {xg) - f'{xg)(Xr - Xg)\ 



(r - ef+' 



where K = ||/'||^ \\x\\}^ and (1-f A)/3 - a > since a < < (1 + A)/3. 

For 9 < < rj introduce the kernel 

G{9, V) J., (C - Of-' iv 0"^' 

ai [a) 1 [2a — 1) 

x^\'"-'(l-g)-" (^1 + (1-9)|3|) dq. (3.6) 



5 



Define for e < a + /S ~ 1 and 9 < < rj < b 



KeAt V) = -Dlr^Dl^-' [G,_ (0, r,)] . (3.7) 



In Lemma [6. 21 we will show that this kernel satisfies 



sup / \K s.t{^, ri)\ d^^dr] < oo. 

is<t<T Js<e<v<t 



0<s<t<T J s<^<ri<t 

Finally, we denote 

A^(x ® y) := / /" Ka.biL ?/)r"-^ {x ® y)^^^ d^drj. (3.8) 

We are ready now to define the integral J^^ f{xr)dyr- 

Definition 3.1 Let {x,y,x ® y) be a (3-Holder continuous multiplicative func- 
tional on M™ (8) W^. Let f : R™— > be a continuously differentiable function 
such that f is X-Holder continuous, where A > ^ — 2. Fix a > and e > 

such that 1 — (i < a < 2(3, a < ^^^^ and e < a + P — 1. Then, for any 
0<a<b<Twe define 

f f{xr)dyr = (-1)" ^ r5^+/.(x)(r)i?il"y^_(r)dr 

m d „i} 

+ EE / Dl%-'d,f,ix)ir)A''^{x^®y^)dr. (3.9) 

Notice that if y is /3-H61der continuous, the fractional derivative Dlz'^yb-ir) 
is well-defined because 

(a-r) 

and (3 + a — 2 > — l.The following theorem asserts that this definition is coherent 
with the classical notion of integral and will allow us to deduce estimates in the 
Holder norm. 

Theorem 3.2 Suppose y : [0, T] M'' is a continuously differentiable function. 
Let X : [0, T] — » R™ be a (3-Holder continuous function and let x ® y be defined 

by {x ® yYg^ = — x'^^ {U'')'^ "^C- Assume that f satisfies the assumptions 

of Definition Wn\ Then, the integral J^^ f{xr)dyr introduced in 1^3. 9\} coincides 

with J2'Li Ja hi^r)iy%dr. 

Proof. To simplify the proof we take m = d = 1. From (|2.8|) and H3.5|l we 

get 



f{xr)y'rdr = {-ir D2+f (x) {r)Dlz''y,^{r)dr 

J a 

= i-ir f' D2+f{x){r)Dlz''y,Ar)dr + A,, 
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where 



Ao = 



r (1 - a) J, J, (r - 6) 



r(i-«) Ja 

So, it suffices to show that 



^^-^Dl-_'^y,.{r)dT\de. (3.10) 



e (r-e) 



A2= DT+-'nx){T)K{x®v)dr. 



(3.11) 



Formula H3.11|l should be first proved for x of class and then extended to 
a general /3-H61der continuous function. Applying H2.3|l . I|2.6|) . and Ij2.5|l we 
obtain 



A. 



a{-lY 



r(i - a) A 



Dl_°'yb-ir) 



f'{xe){xr - xe) 



dOdr 



r(l-a) Ja 



(r - 6')"+i 
f) / 



— Xg 

{r - 6i)"+i 



a(-l)2"- i 

r(i-a) 



^^rv'(x)(e) / 



— 

(r - ey 



where /"'^ denotes the fractional integral applied to a function of 9 and a similar 
notation is used for fractional derivatives. Now 



ja.r x2a — 1 



(r - 6')"+i 

(-l)l-2« rr 



T{2a - 1) 



T{2a - l)r(Q;) Je<e'<^<r'<r 



{r - 6")"+i 

(r - r')''-^{9' - ef'-^ir' - e'^^^dx^dO'dr' . 



Thus 

A2 = 



1 



r(2a-l)r(a) A 

l-b 



Dix-'fixm 



X 



(ry - r') 



9<e'<^<r'<ri 



(3.12) 



(/ - 0')-"-^d0'dr' ) dx^dy^de 
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Making the change of variable = w and using formula 3.196 in Gradshteyn 
and Ryzhik we obtain 

\v ~ r')"~\r' ~ e'y"~^ dr' 



a 



or {v-9'y 



and substituting this expression into l|3.12|l yields 



b 



r{2a - l)r{a)a 



I0<£,<ri<b 

Using 1)3. 6(1 we get 



= G'(6l, ?7)r(2a - l)r(a)a. 
Hence, 

A2 = f Dl%'^f'{x){6) I G{e,^,'n)dx^dy^de. 

Applying the two-dimensional fractional integration by parts formula (j2.9|l 
to ip{^, 77) = G{9, if) and rj) — {x ® y)^,ri and using H3.7|l we obtain 

Je<^<ri<b Je Je 

xT"-%x(giy)^,nd^dr] 

= I f Ke.b{i.v)^''~'{^®y)i..vdidr^- 
Je Je 

Thus 



^2= / Dl%-'r{x){r) f r KrAi,^)r^-^x^y)^^,d^d7jdr. 
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This proves the theorem. ■ 

For any {s,t) G A, and given a /3- Holder continuous multiplicative functional 
{x,y,x ® y), we define 



\\x®y\\s,t,i3 



\Xr " Xg\ 

s<e<r<t \r-9\P 

\{x®y)g^^\ 



™P I a\2B 
<e<r<t \r-e\'^P 



(3.f3) 
(3.14) 



We also set = and = \\x (g) yWo j,^^. Also, \\-\\,^t^^ will 

denote the supremum norm in the interval [s,t\. In the sequel, k will denote a 
constant that may depend on the parameters f3, a. A, e and T. 
The following estimate is useful. 

Proposition 3.3 Under the hypotheses of Definition \S. l\ we have, if b — a < 1 



f{xr)dyr 



< k\ 



a,b,0 



ll2/lla,f,,/3 + fc(ll^®J/|la,fa,/3+ll^llaA^lly|la, 



&,/3 



(ll/'lloo + II/'IIa ll^lla,M (b - (b " ar-'^.{3.15) 



Moreover, if the second derivative f" is X-Holder continuous and bounded, and 
{x,y,x®y) is also a (3 -Holder continuous multiplicative functional onR™®R'', 
then 



f{xr)dyr - I f{xr)dy 

a,b,f3 

< kHi\\x - il|a,6,oo + kH2\\x - X\\a,b,l3 + kHslKx - X)(E) ?/||a,b,/3,(3.16) 

where 

Hi - mup (ll/'lloo + ii/"iIa piu.m {\\x\\ib,p + mib.p)) (b ~ ar'+'^ 

+ (II/"IIoo + ||/"||a(||x||^,,,^ + PII^,,,^) (b-af^) 

X (lla; (8, y\\a,b,p + \\x\\a,b,P \\y\\a,b,p ) (b - af-^', 
H2 = \\f"\\^. \\y\Ufj {\\x\Up + mUp) (b - af('+^'^ 

+ liriL (ll^ ® vLap + MaA, \\y\Up ) (6 - «)'^-'^ 

+ (ll/lloo + ||/||A||i||^,fc,^(6 - a)"') \\y\l ,^^ (b - a)^-2e 
H3 = (ll/lloo + ||/||A||i||^,,,M6- an (6 -a)'^-^^ 



Remark: In H3.15|l we can replace 
ll/'(^)lla.fc,oo' respectively. 



and ll/'lL by \\f{x)\U^ and 
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Proof. First we have, for any r G [a,b] 

|i5:^+/(2;)(r)| <fc(|/(x.)|(r-a)-" + ||/||J|x||^_,,^(r-a)'^-") , (3.17) 



D:+fix)ir)\ < k [\fixr)\ (r - a)- + (r - af+'^''~'^) , (3.18) 

and 

\Dlz'^yt-ir)\<k\\yl.^,^f,ib-rr+^'K (3.19) 
The expression (16. 1() of T yields 

{x ® y)^^,\ < k {\\x ® + \\x\l ,^p (6 - a)2/3+2-2-2e, 

(3.20) 

Consequently, from (|3.8(l Lemma [6.21 we obtain the estimate 

\Ki{x ®y)\<k {\\x ® yll^ + ||x||^ \\y\\a,b,p) {b - 
Thus 



\2/3+2a-2-2e 



(3.21) 



f{Xr)dyr 



<k\ 



X (in ® 2^lla,h,^ + MaAP WyWaAp) " 



-iA/3-2a+l 



2/3+2q-2-2s 



c?r. 



Therefore, we obtain 



f{xr)dy. 



<k\ 



\\y\\a.,,p(P-^Y 



+k Wf'h (MaAP \\y\Ufs + 11^ ® vK.Ab - o^r^^''^- 



2e 



and this implies (|3.15|l easily. 

Note that for any a < 6 < r < b we can write 

f{xr) - f{xe) - f'{xe){xr - xe) - [f{xr) - f{xe) - f\xg){xr - Xg)] 
[f'{xe + z{xr - xe)) - f'{xe)] [xr - xe -Xr + xe)dz 





+ / [f'{.X0 + z{Xr - Xe)) - f'(xe + z{Xr - Xg)) - f'{xg) + f'{xg)] {Xr - Xg 
JQ 

ai - a2. 
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We have 

Kl < \\f"\\^\\x\\aMAU-x\\aM{r-Of^. 

For the term 02 wc make the decomposition 

02 = {Xr -Xe) I I [{Xe - Xe) + Z {Xr -X0-Xr+ xe)] 

Jo Jo 

xf"{xg + z{xr — Xe) + t{x0 — Xe) + tz{xr — Xe — Xr + xe))dtdz 

-{xr -xe){xe - Xe) / f" {xe + t{xe - xe))dt 
Jo 

= {Xr - Xe - Xr + Xe) 

X I I zf"{xe + z{xr — Xe) + t{xe — xe) + tz{xr — xe — Xr + xe))dtdz 
Jo Jo 

-{xr -xe){xe - xe) / [f"{xe + z{xr-xe)+t{xe-xe) + tz{xr-xe- 
Jo Jo 

-f"{xe + t{xe - Xe))] dtdz. 



Thus, 

\a2\ < \\f"\\MaA0\\^-x\\a,,Ar-e)^^ 

+ Iir ||A||^IUA/3lk - 5|U,6,oo {\\x\\i,,0 + \\x - x\\i,^p) (r - e)P^'+>^\ 

As a consequence, 

\f{xr) - f{xe) - f{xe){xr - xe) - [/(ir) - f{xe) - f'{xe){xr - xe)]\ 
< kh{r - e)^^^+^\ (3.22) 



where 



h = Wf'Woo {\\x\\a,b,l3 + \\x\\a,b,^}\\x - X\\a,b,f3 

+ ||/"||API|a,M(l|a^LV + l|iL\/3) \\x- S:\Uc 



On the other hand, we have 



Dl^-'f{x){r)-D',l-'f{x){r) 
1 (f{xr)-rixr) 



r(2-2a) t (r-s)2"-i 

r [nXr)-f'iXr)-f'{xe) + nXe)] ,.\ 

Using the decomposition 

f{Xr) - f{xr) - fixe) + fixe) 

= / f"{Xr + t{Xr - Xr)){Xr - Xr)dt - f" {xe + t{xe - xe)){xe - xe)dt, 
Jo Jo 
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we obtain 



\Dll-^nx){r)-Dll-'nS:)[r)\ 
< k{r - sY-^'-Wf'W^Wx - x\l.r,oo + k{r ~ sf'^''+^\\f"\\^\\x - x\Ur,f^ 

+ k{r - lU + ms,r,p) \\X - i\\s,r,oo 

= fc/2(r-,s)i-2", (3.23) 
where 

h = (iiriioc + II/" 11a (iix||^,,,^ + pii^,,,^) (b - af') \\x - x\u,^ 

+ \\f"\\^\\x-S:\\a.bAb-ci)^- 
Now using l|XT^ . (E!23), we obtain 

l-b 



[f{Xr) - f{Xr)] dyr 



<k\\y\L 



b,0 



X / \f{xr) - fiir)\ {r - a)-" {b - rT+P-^dr 



+kh / (&-r)i-2" \Kl{x®y)\dr 

l-b 



+ / \Di%-'f{S:){r)\\K',{[x-x]®y)\dr. 

J a 

Finally, using H3.21|l we get 



[f{Xr) - f{Xr)] dyr 



<k\ 



lo,b,/3 



(||/'||oo(& - a)'^||x - i|U,^oo + h{b- a)^(2+A)^ 

+ kl2 {\\X ® y\\a,b,0 + \\x\\a.bAy\\a,b,p) (& - af^-^'' 

+[ii/ii.3o+ii/iiAii£iii,,^(&-«n 

X {||(a; -x)®> y||a,6,/3 + ||a; - i||a,6,/3||2/||a,b,/3} (b - a) 



(3.24) 



2/3-2e 



This implies (|XTB|) . ■ 

The following corollary is the direct consequence of the proposition. 

Corollary 3.4 Assume b — a < 1. Under the hypotheses of Definition \3.1\ if 

{x,y,x(®y) is also a (i-Holder continuous multiplicative functional onW"' , 



12 



have 



f{Xr)dyr 



< k 11/11^ \\y - yW^ ,,^^ + fc (||a; (y - y)\\aA0 + MaAP " y\\a,b,p) 
X (ll/'lL, + II/'IL Ikll',.,/, ib ~ - a)^^'^ (3.25) 

On f/ie other hand, if the derivative f" is X-Holder continuous and bounded, 
{x, y, x(S)y) is another (3-Hdlder continuous multiplicative functional on IR™®R'', 
and f is another function satisfying the hypotheses of Definition VS.lX then 



f{xr)dyr - / f{xr)dy.. 



a,b,l3 



< kH(\\x - i||a,fc,oo + kHl\\x ~ x\\a^b,p + kHl\\{x - x) (g) y\\a,b,i3 



+k 



f-f 



x( f-f + f- f ml bpib- ^)'') (b - af-'^. (3.26) 

The estimate ()3.25(l implies that for a fixed x, the mapping (y, a: ?/) — > 
/ f{xr)dyr is continuous with respect to the /3-norm. As a consequence, if y" is 
a sequence of continuously differentiable functions (or Lipschitz functions) such 
that 



\\y-y"\\0 

\x(gy-x®y''\\p 



as n tends to infinity, then 



f{xr)dyr - / f{xr)dy 



0. 



(3.27) 



Hence, the integral J f(xr)dyr introduced in Definition 13 . II does not depend on 
the parameters a and e, and it coincides with the classical integral / f{xr)y'rdr 
when y is continuously differentiable. 

Set fl — — for i = 0, 1, . . . , n. If y is /3-H61der continuous, the sequence of 
functions 



Vt =yo 



1(0} (0 + E W h-i + ^ - ^^-i) - y^-^ 



converge to y in the /3'-norm for any (3' < (3. Assume that the multiplicative 
functional (a^r ~ Xs)dy" converges in the /3'-norm as n tends to infinity to 
{x (8) y)g J. Then 13.2711 holds with /3 = /?'. In particular, this means that 



f{xr)dyr = Jm^X^ ^ i^j^ J{xs)ds^ {ytY - yti_^)- 



(3.28) 
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For any p > 1, the p variation of a function x : [0, T] ^ R is defined as 
Ye.rp{x) = Slip ( ^ \x{t^) - xit^)]"] , 

where n — {0 — to < ■ ■ ■ < tn — T} runs over all partitions of [0,T]. Notice 
that 

Yavi/pix) < \\x\\fj. 

Then, for any /3-H61der continuous multiplicative functional {x, y, x®y) on M™(8) 
M'^ and any function / satisfying the hypotheses of Definition 13.11 the integral 
/o f{xr)dyr coincides with the integral defined using the ^-variation norm (see 

This implies that jj" f(xs)dya is given by the limit of the Riemann sums 
of the form 

Jo KHo^ 

where n = {0 — to < ■ ■ ■ < tn — T} runs over all partitions of [0, T]. 

In order to handle differential equations we need to introduce the tensor 
product of two multiplicative functionals: 

Definition 3.5 Suppose that {x,y,x (8" y) and {y,z,y (g) z) are (3-Holder con- 
tinuous real valued multiplicative functionals. Then, for all a < h < c, we 
define 



{x®{y(®z)\ = Aa,r,bA^^z)Dl_°'yt-{r)dr 

^ ' / a.b ,l„ 



r(2 — 2a) Ja<r<^<ri<b 

r"^" (a: y)^ ,^ (z, - Zr) ~ {xr ~ Xg) F"-" {y (g> z)^^^ 
^ ' (r - a)2"-i 

where 



drd^dr/, 



(-1) I [Xr - Xa){Zc - Zr) , f (z^ - Zg) (xg - X^) 

Aa,r,b,c(x,z) ^ — +a -—r^ d9 

ri-a \ r-a)« (^r - 9)^ 



We have the following result. 

Proposition 3.6 // the function y is continuously differentiahle and for all 
a<b 

iy ® ^)aM = i^b- Zr) y'rdr, 

J a 

i^^y)ab = / {Xr - Xa)y'rdr, 
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then 



(x (g) (y (g) z). ^) = / {Xr - Xa){Zc - Zr)yl. 
\ ' / a,b . „ 



dr. 



Proof. We are going to use formula H3.9|l with m = 2, d — 1, f{x, z) — xz 
and the functions xt — Xa and Zc — zt . In this way we obtain 



{xe - Xa){zc - ze)dye 

(-1)" I {Xr - Xa)(Zc- Zr) f {Zr - Zg){xg - Xr) 



xDlz''yb-{r)dr 
1 rb 



(r - a)" 



r(2 - 2a) J, 



+ [2a - 1 / — 

[r - a)^" ^ Ja {f- 0) 



{r-ef 



Zr - Zg 



\2a 



dO 



dO 



1 



r(2 - 2a) A 



X.,f,(e,ry)r"-^ {x®y)^^^d^drjdr 
b r z*^ 

(r-a)2"-i (r-0)2" 



— Xg 



de 



X.,6(^,77)r"-^ {y®z).d^dridr. 



and this completes the proof. ■ 

It is easy to obtain the following estimate 

Proposition 3.7 Suppose that {x,y,x®y) and {y,z,y(!S)z) are (3-Hdlder con- 
tinuous real valued multiplicative Junctionals. Then, for any a < b < c we 
have 

(x (y O ^).,c)^ ^ < k (||yL,h,/3 Ma,b.P Ma.b,(3 

+ fc \\4a,c,f3 {\\y\\a,b,0 Ma,b,P + WV ® ^\\a,b,f3 ) (& " " a)^ . 



If 6 = c we write (x(8)(y(g)z)jj) =(x(8)j/(8) z)^ j^. If the functions x, y 

V ' / a,b ' 

and z are continuously differentiablc, then 



(x (g) y (g) = / x'^y'gz'^drdOda. 

J a<r<e<(7<b 

Define 

II ^ ^ II |(x(8)t/®^)g I 
||a;®y®z||^^^= sup , .,3^ ■ 
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Then, Proposition 13 . 71 implies that 



{\\v\\a,b,p \Ma,b,P Ma,b,0 + Ma,b,p WV ® Aa^ 



l|a^lla,fc,;3b®^llaA/3 



(3.29) 



Proposition l3 . 7l also imphes that {x, {y ® z). ^ , ^a; O (y z), ) is a /3-H61der 
continuous functional on the interval [0, c]. As a consequence, if / satisfies the 
assumptions of Definition 13.11 we can define the integral J^^ f{xr)dr {y ® z)^ ^, 
for all a < 6 < c. The following estimate for this integral will be needed to solve 
differential equations. 

Proposition 3.8 Suppose that {x, y, x®y) and {y, z, y®z) are (3-Hdlder contin- 
uous multiplicative junctionals on (8)R''. Let f : R'' be a continuously 
differentiable function such that /' is \-Hdlder continuous, where A > ^ — 2. 



Fix a > and e > such that 1 — f3 < a < 2(3, a < 
Then the following estimate holds 



A/3+1 



and e < a + P — 1. 



1 



f{Xr)dr{y » z)r,r, 



< k 



An 



+Ba,b\\x®y\\^^^^0 {b~a) 



/3-2£ 



where 

Aa,b 

and 



\\y z\\a,b.P + \\y\\a,bAAaAfj) 



(3.30) 



(3.31) 



ll/lloo + {MaAP ll/'lloo + II/'IIa MIa, (b - a)^f') {b - af-'^ 



Ba.b = \\z\Up (II/'IU + II/'IIa \\<b.p (b - a)'^) ■ (3.32) 

Proof. To simplify the proof we will assume d = m = 1. From (|3.1() it is 

easy to see that 

l^"" ® y^-4aA0 - (""^ ® ^""■'■^ "''"''■^■^ WyWaAp) - (3.33) 

and from Proposition 13 . 71 we have 

llx ® (y Z), f, \\a,b,f3 < fc( |NIL,6,/3 \\y\\aA0 ll^lla,fc,/3 

+ II^La^^ \\y ® 4a,b,^ + \\4a,b,^ Ik ® 2/IL,M ) - • (3-34) 



16 



From (26), 1(223, and we obtain 




f{Xr)dr{y z)r,t 



< k\\f\\oo\\{y^z).4a,bAb-a)'^ 

+k[{\\x (y z).^b\\a,b,i3 + ||a;||Q,b,/5||(y z).^b\\a.b,f3) 
X (ll/'lloo + ||/'IU||x||^^,,^(6 - a)^0)]{b - af^-^^ 

< fc||/l|oo (112/ 2||aA/5 + l|yL,f,,/3lNLA/3) - 

+ fc (l|2;|la.b,/3 ll2/lla.b,/3 \\4aAI3 + Ma^b.p WV 4a,b,l3 + MaAP ® yL^b.p) 

X (ll/'lloo + ||/'IU||2;||t,,^(6 - af^)]{b - a)3^-2^ 
which implies the desired result. ■ 



4 Differential Equations Driven by Rough Paths 

Let y : [0, 1] ^ R'* be a /3-H61der continuous function. Suppose that {y\y^ ,y^(S> 
y^) is a /3-H61der continuous multiplicative function, for each z, j = 1, . . . , d. We 
aim to solve the differential equation 

Xt=XQ+ I f{xr)dyr, (4.1) 
Jo 

where / = M™ ^ M™'^. 

Formula 1)3. 9|) and Definition 13.51 allow us to transform this equation into the 
following system of integral equations: 

Xt = XQ + {-ir f D^+f (x) {s)Dl~''yt^{s)ds (4.2) 
+ Dll-'f {x) (s) ^* J' KoAt '?)r"-^ (x ® y\^^ d^dr,ds, 



t 

l-a 



i^^y)s,t = (-1)" / i??+/(x)(r)A^-"(y®2/)..,_(r)dr 



Dlr'f i^)ir) (4.3) 



K 



^r (C, ?/)r" ' (x®{y®y).^ d^diqdr 



Theorem 4.1 Let y : [0, 1] ^ K*^ he a (3-Hdlder continuous function. Suppose 
that (y', y-' ,y^ ®y-') is a real valued f3-Hdlder continuous multiplicative function, 
for each i,j = l,...,d. Let f : R™ M™'^ be a continuously differentiable 
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function such that f is \-Hdlder continuous, where A > ^ — 2, and f and f 
are bounded. Set 

Pf :=II/I|oo + ||/'||oo + ||/'||a. 

Then there is a solution to Equations j^-i^ - such that {x,y,x y) is 

a P-Holder continuous multiplicative functional. Moreover, for any 7 > ^ the 
function x satisfies the estimate 



sup |a;t| < IxqI + T<2kpf 

0<t<T 



\y(Sy\\i} 
\\yh 



V 1 



(4.4) 



where k is a universal constant depending only on [3 and 7. 



Proof. To simplify the proof we will assume d = m = I. The proof will be 
done in several steps. 

Step 1. Fix a > and e > such that 1 - /? < a < 2/3, a < 
e < a + /3- 1, £ < f , and (1 - 2e)/(/3 - 2e) < 7. 

We will write the Equations (|4.2|) and (|4.3|) in the compact form 

X = <^i{x,y,x®y), 
x(g)y = <^2{x,y,v®y,x®y). 

Consider the mapping J : {x,x ®y) (Jix, J2 {x ® y)) defined by 



JlX 

J2 (x y) 



<Pi{x,y,x(g)y), 
^2{x,y,y®y,x®y). 



We need some a priori estimates of the Holder norms of Jix and J2 {x ® y) in 
terms of the Holder norms of x and x y. From 13.15|l it follows that 



< k 



m\\oo\\y\Up + {\\^®y\\s,t.p + Ms,uf>\\y\Up) 



X (ll/'lL + II/'IIa Mit,p it - ^)''') it - sr-'']-i^-5) 

On the other hand, Proposition 13.81 implies that 

II J2 (x ® y)\\,^t^p < k \As,t + B,,t \\x ® y\\^ ^ p {t - sf-^' 



(4.6) 



where Ag^t and Bg^t are defined by H3.31|l and H3.32|) . respectively. 
Step 2. Set 



a{y) := 2kpf 



\yh + 



\\y®v\\(} 
hh 



V 1 



l/(/3-2£) 



where k is the constant appearing in formulas 14. 5|) and 1)4. 6|l . Suppose that 

1 



< t - s < 



aiy)' 



(4.7) 
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Then, the inequahties 



Ms,t,p < 2kpf\\y\y (4.8) 
{t~sf-'^\\xC=^y\\^^,^^ < \\yy (4.9) 



imply that 



\\Jix\\s,t,p < 2kpf\\yy (4.10) 
{t~sr-^^\\Mx®y)\\,^,^p < \\yy (4.11) 

In fact, from the definition of a{y) and 14.8|l we deduce 

{t - sf-^^\x\\s,t,p < I. (4.12) 

By the definition of Bg.t and As^t we have 

Bs,t < (ll/'lloo + ||/'IU)||2/||/3 < PfWvh < ^^''^fc'^'''^ (4.13) 

and 

As,t < (||/||oo + ||/'||oo + ||/'||A)(||y®z/||^ + ||y||; 

< hh- (4.14) 

Therefore, substituting 14.13|l and (|4.14() into (|4.6|l we obtain (|4.11() . Finally, 
from we get (QUI) . 

Step 3. We can now proceed with the proof of the existence. Let be a 
natural number such that = S < ^j^^/j. We partition the interval [0,T] in A'' 
subintervals of the same length and set ti = i ~ 0, l,...,iV — 1. We will 
make use of the notation \\x\\^ — ||a;||tj_i,ti,/3, and \\x (X> = ||a;(8)2/||ti_i,ti,/3, for 
i = 1, . . . , N — 1. From Step 2 we know that if that x and x y satisfy 

||x||. < 2kpf\\yy 

for any i = l,...,iV — 1, then the same inequalities hold for Jx and Jx y, 
that is 

||Jia;||, < 2kpf\\yy 
]\J2{x®y)l < \\y\\f,S-^''-'^l 

Consequently, there is a constant Ci such that 

\\J^xy + \\J^{x<E>y)\\p<Ci. 

This implies that the sequence of functions Jf x is equicontinuous and bounded 
in C'^. Therefore, there exists a subsequence which converges in the /3'-H61der 
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norm if /?' < f3. In the same way, there is a subsequence of J^(a; y) which 
converges in the /3'-H61der norm. The limit {x,x'S)y) defines a /3-Holder contin- 
uous multiplicative functional (x, y,x®y). Using the continuity of the solution 
in this norm it is not difficult to show that the limit is a solution. This implies 
the existence of a solution, which satisfies H4.8|) and H4.9|l . 

Step 4. Let us now prove the estimate (|4.4ll . By step 2, the solution we 
have constructed satisfies the estimates (|4.8|l and (I4.9|l if H4.7|l holds. Then it 
follows that for any r G [s, i] 

sup \xr\ <\xs\ + {t~- s)''||a;||,,t,;3 < + (t - sf. 

r£[s,i\ 

Since the interval [0,T] can be divided into [T/r] intervals of length r = -^^^ 
the inequality H4.4|l follows. ■ 

Theorem 4.2 Let y : [0, 1] he a /3-Hdlder continuous function. Suppose 

that (y', y^ ,y^ ®y-') is a real valued (3-Hdlder continuous multiplicative function, 
for each i,j — l,...,d. Let f : M."^ — > R""' be a twice continuously differentiahle 
function such that f" is X-Holder continuous, where X > ^ — 2, and f , f and 
f" are bounded. Then there is a unique solution to Equations j^.i^ - such 
that (x,y,x (®y) is a [3-H6lder continuous multiplicative functional. 

Moreover, if x satisfies xt = xq + f{xr)dyr and y verifies the same hy- 
potheses as y, then 

sup \xt - £t\ < C {\xo - io\ + \\y - yy + \\y (g> (y - y)y} ^ (4-15) 

0<t<T 

where C depends on \\y\\/3, \\y®y\\p, (3, \, and pf, and where 

pf - ii/iloo + ii/iu + ii/'iloo + ii/'iu + II/" Hoc + iinu ■ 

Proof. To simplify the proof we will assume d = m = \. Notice that 
uniqueness follows from the estimate (|4.15|l . So it suffices to show this inequality. 
We fix s < i such that t — s < , where /3 is defined as follows 



Piy) = 2kpf 



|y||/3 + ll2/ll| + lly®y||/3- 



\\yh 



l/(/3A) 

VI) . (4.16) 



The constant k appearing in the definition of /3 will be chosen later. We choose 
a > and e > such that 1- P < a < 2P, a < and e < a-\- fi - 1, e < ^. 

We also assume that the solutions x and x satisfy the following inequalities: 

{t-sr-'^\\xy<l, (4.17) 

{t-s)^-'^\\x\\p<l, (4.18) 
it-sf-^^\\x®y\\^<\\yh. (4.19) 
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Our first purpose is to estimate the Holder norm \\x — x\\^ f p. We can write 

f{xs)d{ys - Vs) 



\X-X\\s,t,0 < 



s,t,l3 



[f{xs) - f{xs)] dy. 

— h,s,t + l2,s,t- 

The term Ii^s,t can be estimated using and we obtain 

h,s,t < k[Hi\\x - i|U,t,oo + H2\\x - S:\\sj,p + HsWix - x) <E) y\\s,t,fi]^ (4-20) 
where 

Hi - l|ylL^,,^ (ll/'lloo + liril, PIU,,,^ + \\x\\lt,p)) {t - sf^'^'^ 

+(ii/"iioo+ii/"iu(iixii^,,,^+Pii^,,,^) (t-^n 



Ho = 



\\r\L\\y\Up(Ms,t.p + ms,t.p){t~sf('+^^ 

+ iiriloo (ii^ ® yh^fi + MsA^ \\y\\s,t.,, ) it - ^r^-'^ 

loo + WfWxMitAt-^)"''' 



s,t,l3 



it~s) 



/3-2s 



and 



H3 = (ii/iloo + wfhmitAt - ^)'^) it - ^f-'^- 

Then, using the inequahties H4.17|l . (|4.18|l . and H4.19|l we get the following esti- 
mates 



(4.21) 
(4.22) 
(4.23) 

It remains to handle the term \\{x — i) ® y\\s,t,f3 in H4.20|l . To get estimates for 
this term we apply again the inequality H3.16|l and we have 



Hi < (ll/'lloo + 2||/"||oo + 6 liriU), 

H2 < 112/11^3 (3 ||riL+ ||/||oo + ||/||A)(t-s)''\ 

H3 < (II/IIoo + II/IIa). 



\{{x - i) (g)y)sA 



IfiXr) - fiXr)] driy ® y)r.t 



< k{t-sf Hi\\x~S:\\s,t.,oo+H2\\x-x\\s,t.,p 



+H3\\{x - i) (g) (y (g) y).,t\\s,t,0\ , (4.24) 

where 

Hi = ||(y ® y).As,t..0 (ll/'ll- + liriU W^lUp {Ml,, + Mlj,)) it sr^^'^'^ 
+ (||/"||oo + ||/"||a(||x||^,,,^ + PI|^,,,^) {t-sfl 

X (llx ® (2/ ® y).,t\\s,t,p + Ms,t,p Wiy ® y);t\\s.t,p ) it - 
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H2 = liriL Wiy ® y);t\\s,t,p + Ms,tj3) {t - sf^'+^^ 

+ II/" lloo (ll^ ^ (y ^ y)-As,t,0 + Ms,t,f> Wiy ® y).,tlL,,,^ ) (t - sf^-^^ 
+ (ll/IU + WfUmitdt - 'f^) W^y ® y)-As,t,, - ^f'^'^ 

and _ 

Using lESH, (EIZI), and (j^l^ we get the following estimates 



^1 = (lly®ylL,t,0+ \\y\\s,t,p^ 

X (ii/'iloo + iiriu piim,/3 (ii^iIm,^ + ii^iIm,^)) (t - 
+fc (iiriioo + II/" iia (ii^ii^,,,^ + pii^,,,^) (t - sf^) 



X (lkL,t/3 \\y\\lt,0 + Ms,t,fi \\y » y\\s,t,p + \\y\\s,t,0 11^ ® ylL,*,/?) - ^) 



< k 



{\\y ^ vh + Ml) (ll/'lloo + ll/"lloo + liriL) {t - sf, 



20-2e 

(4.25) 



and 



+fc|l/"lloo(ll2/llM,/3ll^®y|lM,/3 (4-26) 
+ \\A\s.t,0\\y®v\\s,t,0 + \\x\\s,t,p hWltji ){t-sf^-^' 

+k (ll/lloo + ll/IUI|i|lL,/5(i - ^)^^) {\\y ® y\U, + Mlt.p) (t - ^r^-'' 
< k (ii/iu + ll/IU + lir lloo) (lly ®yh+ \\y\\i) it - ^f- (4-27) 

On the other hand, from (|3.34() we get 

||(x - x) ® (y y).,t\\s,t,/3 < fc( llx - i||^ \\y\\l 

+ Ik - s:\\s,t,p \\y ® yll/3 + \\y\\p ll(^ - ® yh.t.p ) - ^)''- (4-28) 

Thus, substituting iITTtIi . and into (jO^ yields 

||(a; - i) ® yl ^ ,, < k{t - sf [ (||y ® + \\y\\)^ 
x(||/'||oo + ||/"||oo + |iriU)||x-i|U,*,„. 

+ (ll/lloo + ll/IU + lirilj (||y ® y\\p + ||y||^) Ik - 

+ (ll/lloo + ll/IU) 

X (ik - ^IL,t/3 llyll^ + Ik - ^IL,t,/3 lly ylU + hWp -x)® yW^^^^p 
k{t - sfpf {\\y ® y\\i3 + llyll^) (Ik - ^lls,t:oo + Ik - ^lU,*,/?) 



< 



+k{t - sf (ll/IU + ll/IU) II2/IU ll(a; - i) ® 2/IU,,,^ . 
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The condition t — s< if the constant in P{y) is chosen in an appropriate 

way, impHes that 



fc(i-^^)''(||/||cc^ + ||/||A)||2;||^< 



Hence, 



Wi^ - i) «) y\\,^t^p < k{t-s) ^pf {\\y ® y\\f, + (||.t - + \\x - i\l,t,0) 

(4.29) 

Substituting (jO^ . llOT|l . (jO^ and lH:^ into (fi:^ yields 

-fi,s,t < fc/5/[||y||^ ||a; - .Tll^^t^oo + II2/II/3 - i||s,t,/3(i - s)^^ 

+ (\\y ®y\\D+ bll^) (11-^ - i\\s,t,oo + - i\\s,t,0) [t - s)'']. 

Again, condition t — s < 1/ (3{y), if the constant in (3{y) is chosen in an appro- 
priate way, imphes that 



(4.30) 

For the term l2,s,t we have the following estimates, using (|3.25|l 

l2,s,t < k ll/IL ||y - y\\s^t..p + k (ik ® (y - y)\\s,t,0 + Ms.t.,p h ~ ylL.t./j) 
x{\\f'\L + \\f'hMs,.p it-^'^)'') [t-sf-^^ 

< kpf \\y + k\\x®{y- y)\\s,Ui (II/'IL + II/'IIa) " ^f-^^^M) 

In order to estimate \\x ® [y — y)\\g f p we make use of Proposition 13.81 and we 
obtain 



\\x®{y-y)\\s,t,p 



1 



f{xr)dr{y ® (2/ - y))r,n 



< k 



A,^t + Bs,t\\x(E)y\\^^^^U {t-sf-^' 



(4.32) 



where 

As,t 



= {h ®{y- y)\\s^,p + llylL,t,/3ll(y - 2/)IL,t,/3) 

X il/lloc + (ll^lL,,, II/'IL + II/'IIa II^IIm'^* ~ ^)'^) - '^)''"'^ 
< Pf {\\y ®{y- y)h.t,fs + \\y\\s^t,My - y)!^,*,/?) > 



(4.33) 



and 



Bs,t = 11(2/ -y)IL.,,(ll/'lloo + II/'IIa II^IIm,/3(*-^) 
< ll(2/-y)llM./3(ll/'lloo^ 



A/3 



(4.34) 
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Substituting and (jOIjl into yields 

\\x^iy-y)\\s,t,i3< kpf[\\y®{y--y)h + \\yh\\iy"y)h]- (4-35) 
Finally, from (|4.35() and H4.38|l we obtain 

l2,s,t < kpf \\y -y\\p + kp) [\\y ® {y ~ y)\\0 + Ms.tMv ~ Mfi] (t ~ ^f'^' 
< kpf\\y-y\\p + kp}\\y(g>(y-y)\\p(t~s)^'''- (4.36) 
Now from (|4.30|) and H4.36|l we get 

\\^-^\\s,t,f} < kpf\\y\\f^\\x~S;\\,^t,^ + ^\\x-x\\sA,i3 

+kpf \\y - y\\^ + kp}\\y ® (y - y)y{t - sf-^'. 

Or 

+kpf \\y-y\\(} + kp)\\y ®{y- y)\\p{t " s)''^'^(4.37) 

Notice that 

\\x - S:\\s,t,oa <\xs ~ S:s\ + {t - sY\\x - x||s,t,/3- (4.38) 

Hence, 

\\x-x\\s^t^p < kpf\\y\\p[\^s-Xs\ + {t-sf\\x-i\\s,t,p] 

+kpf \\y - yWp + kp}\\y ®{y- y)\\p{t " sf'^"- 

And consequently, 

\\x-x\\,tf3 < ^Pf WvWp - S:s\ + kpf \\y - y\\p + kp)\\y® {y - y)\\f3{t - sf-"^" . 

(4.39) 

Substituting into yields 



s,t,oc ^ \Xs - xs\-t{t~ sY{kpf \\y\\p \Xs - Xs\ 

+kpf \\y - y\\p + kpjWy (y - y)hit - s)"-'') ■ (4.40) 

Suppose that y = y. Then, Equation (|4.4(J|) implies that x = -r in a small 
interval [0,6], and by a recursive argument, the uniqueness follows. 
Denote k = and i„ — nn. Set 

Zn = sup \Xs - Xs\ 
0<s<t„ 

Then inequality (|4.4U|) states that 

Zn+i < (1 + kpfKf')Zn + kpfK^Wy - y^ + kp}^^-^'\\y ® {y - y)\\p 
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Therefore 



T/k 

+fc^(l + kpfnPy [pfKf'Wy - ijWp + pjn'^-'^y ® (y - y)h] 



1=0 



This impUes the desired estimate. ■ 

The following corollary is direct consequence of H4.H7|I and (|4.15(l . 

Corollary 4.3 /// is twice continuously differentiable and f" is Lipschitz con- 
tinuous and if x and x satisfy 

Xt = xo+ / f{xs)dys and it = xq + / f{xs)dys , 
Jo Jo 

then 

\\xt - itWp < C {\xo - io\ + \\y -y\\f} + \\y ® (2/ - y)h} , (4-41) 

where we use the notation of Theorem 4-1- 

5 Stochastic Differential Equations 

Suppose that B = {Bt ^ {Bl,Bf, . . . , Bf)} is a d-dimensional Brownian mo- 
tion. Fix a time interval [0,r]. Define 

{B(g,B)^ ^ = [ {Br -Bs) do Br, 
J s 

where the stochastic integral is a Stratonovich integral. That is, 

\{Bl-Blf if i^j 
^^""^^-^'X jliBl-B^dBl if ' 

where the stochastic integral is an Ito integral. It is not difficult to show that we 
can choose a version of {B ® B)^ ^ in such a way that [B, B,B ®B) constitutes 
a /3-Holder continuous multiplicative functional, for a fixed (3 G (1/3, 1/2). 

As a first application of Theorem l3.2l and (|3.28(l we deduce that the Stratonovich 
stochastic integral f{Br)d o Br has the following path-wise expression 

r f{Br)odBr = {^ir f ^0+hiB)r{D}f_^B'T-)rdT 

Jo Jo 

+ EE/ Dl^-'dJ,{B)^Aj{B^C^jT^)dr. (5.1) 
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We can apply Theorem 14.11 and deduce the existence of a solution for the 
stochastic differential equation in M™ 

Xt^XQ+ f f{Xs)dBs, (5.2) 
Jo 

where the initial condition Xq is an arbitrary random variable, and the function 
/ : R™ — > R""' is a continuously differentiable function such that /' is A-H61dcr 
continuous, where A > ;g — 2, and / and /' are bounded. By Theorem 14.21 
the solution is unique if / is twice continuously differentiable with bounded 
derivatives and /" is A-Holder continuous, where A > ^ ~ 2. The stochastic 
integral here is a path-wise integral which depends on 13 and B ® B. 

We have also the stability type results (|4.15() and (|4.41() . In particular, if B"^ 
is a piece-wise smooth approximation of B such that 

\\B'-B\\p and \\B {B^- - B)\\p 

converge to zero with a certain rate, then according to Corollarv l4.3l — X"^!!^ 
will also converge to with the same rate, where 

Xt ^Xo+ f f{Xl)dBl 
Jo 

In particular, this implies that the stochastic process solution of H5.2|l coincides 
with the solution of the Stratonovich stochastic differential equation 

Xt = Xo+ f f{Xs)d°Bf. (5.3) 
Jo 

In this section we will apply these results in order to obtain the almost sure rate 
of convergence of the Wong-Zakai approximation to the stochastic differential 
equation (|5.2|l . That is, we will consider the rate of convergence in Holder norm 
when we approximate the Brownian motion by a polygonal line. 

In order to get a precise rate for these approximations we will make use of 
the following exact modulus of continuity of the Brownian motion. There exists 
a random variable G such that almost surely for any s, f G [0, T] we have 

\Bt - Bs\ < G\t - s\^/^^log{\t - s\-^) (5.4) 

Let TT = {0 = to < ti < ■ ■ ■ < tn — T} be the uniform partition of the in- 
terval [0,r]. That is tfe fc = 0, . . . ,n. We denote by B'^ the polygonal 
approximation of the Brownian motion defined by 

k=0 

We have the following result 
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Lemma 5.1 There exist a random variable Ct,p such that 



(5.5) 
(5.6) 



Proof. Fix < s < t < T and assume that s G [ti,ti^i] and t e [tk,tk+i]. 
Let us first estimate 



/ii(s,t) 



1 



{t-sr 



\B^ -Bt- (BJ - B, 



If t-s > -, then using (|5.4|) we obtain 

\hi{s,t)\ < T-PnP\Bt,-Bt + ^(t-tu){Bt,^,-Bt,) 



Bt 



Bs + - {s - ti) {Bt,^, - Bt,] 



< 4Gr-'5+i/2n-i/2+/3 0og(n/r). 

Ift— s< ^, then there are two cases. Suppose first that s, i € [tfe,</c+i]. In this 
case, if n is large enough {n > Te^/'-'^"^''^) we obtain using H5.4|l 

\h t A\ <, \Bt^B,\ , n \Bt,^,-BtJ ^^ 
\h^{s,t)\ < -^-—^ + -^-—^it-s) 

< G\t - s|3-/3^1og \t - + GT-^/^^login/T) n^~^/^{t - s)^-^ 

< Gr-'='+i/2„-i/2+/3 0og(n/r). 

On the other hand, if s £ [tfc_i,tfc] and t £ [tk,tk+i] we have, again if n is large 
enough 



\hi{s,t)\ < 



< 



< 



1 






— < 


Bt^, - 




1 r 








G 


it 


~s)P 



Btk ~ Bt + — {t~ tk) {Bt^^^ - Bt^) 



\Bt-Bs\ + -{t-s){\Bt,-Bt,^A + \Bt,^.-BtJ) 



|i~s|i/Viog|t-s|-i + 2(t-s)(^) v/iog("/r) 



This proves H5.5|l . 

Now we turn to the estimate of the term 



h2{s,t) 



1 



{t - 5)2/3 



- Bl)dB^^^ ~ J\b: - Bl)dBi 
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for i ^ j (the case i — j is obvious from (|5.5|l . We claim that the there exists 
a random variable Z such that, almost surely, for all s,t 6 [0, T] we have 



fiBl-BDdBi 

J S 



< Z\t-.s\log\t-s\ 



(5.7) 



In fact, it suffices to show this inequality almost surely for all s and t rational 
numbers. If we fix s, the process {Mt,t g [s,T]} 



Mt^ j\Bl~Bl)dBi 



is a continuous martingale and it can be represented as a time-changed Brownian 
motion: 

As a consequence, applying 15. 4|) there exists a random variable Gi such that 

1/2 



\Mt 



Wj^'(_Bj_Bj)2 



< G (^^ (B; - Bl fdu^ y log (^^ {Bl, - 



and again (|5.4(l , applied to Bl^ — Bl, yields 



\Mt\ < G1G2 



(u — s)log\u — s\ ^dii^ ylog^G| y (w — s) log |m — s| ^dii^ 



for some random variable G2 . We have for | i — s | < 1 

y (u - s) log |w - s\-^du ^\t- Q + ^ log \t - s\ 

and this implies easily the estimate (I5.7|l . 
Suppose first that t — s > — . Then 



ft.2(i, s) 



< 



1 


it 






1 


it 










1 



(i - s)2/3 
1 



l\Bl^~Bl)dB^-~ l\Bl~Bl)dBl 

{Bl - Bl){Br - Bi'^ fiBi;- - Bi^dBl 

J S 



{Bl-BD{Br ^Bl-^-Bi-Bi) 



{t - s)2/3 

Ai +A2. 



f [B^ - Bi] dBl 
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Using (|5.4|l and (|5.5I) the term Ai can be estimated as follows 

Ai < G\t - s|i/2-'3 0og|t - - B^Wp 

< CT,/3n'3"^/Vlog"- (5-8) 
For the term A2 we proceed as in the proof of the estimate (|5.7|l . We have 

where is a Brownian motion. As a consequence, using that 

lis - B^Woo < CT,pn-^'^^/\og{n/T) 
(this estimate is proved as H5.5|l ) we get 

1/2 ' 



Ao < 



G 



{t - s)2/3 



^ - Bifdu^ ^ log ^ - 



Suppose now that t — s < We make the decomposition 
1 



(5.9) 



/l2(s,t) < 



(t - s)2/3 
i?l + i?2 ■ 



- + ^V; - Bi)dB- 



Then (|^ yields 

B2 <Z\t-s\^-'^'^\og\t-s\-^ <CT,pn'^'^~^\ogn. (5.10) 
In order to handle the term Bi, assume first that s,i G [tk,tk+i]. Then 

J\bI - B^dBir - ^ (bI,, ~ bQ J\bI ~ Bl)du 
and we obtain 

Bi < CT,p{t - sy~^l^\ogn < CT,pn^'^'^ \ogn. 
Finally, if s £ [ifc_i,tfc] and t G [tk^tk+i] we have 



K Bl_^) f\Bl - Bl)du + [bI^^ bQ fiBl - Bl)du 



Bi < (t^s)-^'^^ 
< CT,0n^~^'^ log 



The proof is now complete. ■ 

As a consequence, we can establish the following result. 
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Theorem 5.2 Let f : R"* — > R"*'' be continuously differentiable with bounded 
derivative up to forth order and let X satisfy 

Xt^Xo+ f f{Xs)dB, 
Jo 

If satisfies the following ordinary differential equation 

X^ - Xo + / f{X^)dB^ , 
Jo 

then for any f3 G (1/3, 1/2), there is a random constant Ct.p £ (0, oo) such that 



\\X~X^\\p<CT^pnP-^'^.J\^. (5.11) 

Proof. The result is a straightforward consequence of Lemma I^TI and The- 
orem ■ 

6 Appendix 

Proof of Lemma l2.2L The fractional integration by parts formula H2.8|l yields 



= (-1) 



b pb 

1-a 



The operators D^^"'^ and ^ commute, as it follows from the following com- 
putations: 



^ ^ [(^ - 0"" V(e, v)] + (« - 1)('7 - 0""'^(C, v) 



r(a) [ 



1-a 



(ry-0^ 

r(a) 04 I h {-n- v ) 



-i?^i;-V(e,77). 
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Thus 



Hence, applying again H2.8|l we obtain (|2.9|) with r°'i/j{£^,ri) given by H2.10|l . 
Remarks: 

1 . Formula (|2.9(l holds if ip is of class in a < ^ < 77 < 6 and 



d^dr] < 00. 



2. Under the conditions of the above lemma, we also have T"ip{£^, rj) 

3. The operator F" can also be expressed as follows. 



T{af 1 



+ (!-«) / (e'-O 



iv - v') 

(77-er-V(e,?7)-(^-?'r~w,r?) 



+ (l-a) 

Tiaf 
+ (!-«) 



« (77 - r/)' 



-drj' 



dC 



iv - v'Y 



drj 



Exchanging the integration order, we see the last double integral equals to 



iv - v') 



-d^'di 
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This leads to the following expression for 

r{a) 

^^'-"^l 

Consider the kernel Ks.t{£,,'ri) defined in H3.7|l . that is, 
where 







(6.2) 



1 



JQ Jo 

(6.3) 

and Cq is given as the coefficient in (I3.6|) . 

Lemma 6.1 Lei 1/2 < a < 1. T/ie function 0(z) defined in 16. '-ip satisfies 
(/>(0) < oo, (/) is decreases to zero as z tends to infinity. If (3 < I ~ a, then 

(j){z) < cz-^ (6.4) 

Moreover, if [3 < 2 — a, 

|0'(z)| < cz-f^. (6.5) 



Lemma 6.2 The kernel Ks^t{S.,v) satisfies 

sup / \Ks.ti^,Tj)\d^dr] <(x,. (6.6) 

0<s<t<T Js<5<,)<t 

Proof. To simplify the notation we omit the dependence on the variable s 
in G(s,^, 77). Also, c will denote a generic constant depending on a and e. We 
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have 



^""^ I (t-.)- "'^"-^V, (.'-,)"-^ '^'^ 
V(t-77)"-^(^-s)"-^ (77' - 77)"-=+! 

7)-G(e,V)-G(e',^ 



Set 



1 /■*G(e,??)-G(^,r;'), , 



A ^ Gii,v)-Gitt)-GiC,v)+GiC,t) 



It suffices to show that 

sup / \Ai\d^dr] < 00 (6.7) 

0<s<t<T J s<i<7j<t 

for i 1,2,3,4. 

Step 1 Suppose i — I. Using the fact that the function cj) is bounded we 
obtain 

Hence, 

\Gitv)\ + m,t)\ 



< c(e-sr^(ry-e^"'(^~ry)-"^ 



and (|5.7h holds for i = 1. 

Step 2 Suppose i = 2. We have 



G{^,v)-G{U)^ -Tri^,y)dy, 



V 



dy 
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and 



-(e-.sr-i(7?-e)"-V' 



where 



Notice that, by Lemma l6. II the function x(^) is uniformly bounded. Hence, 

\A2\ < c(C - sy-' f r {y - 0"-'(ry' - v)-"+'-'dydrj' 

< c(e - sf-^r^ - f r (r/ - yf^-\y - r,)i-Hydi 

J q J r] 

which imphes that (|6.7I ) holds for i — 2. 
Step 3 Suppose i = 3. We have 



G(e,^)-G(^',77) = J -^{x,ri)dx, 



and 



-(a-l)(^-,sr-i(r7-a"-' 



where 



= (e-sr-^(r,-er-'(a-i)0 

7(z) = (l-a)(/)(z) + (l + z)(/.'(z). 
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By Lemma l6. II the function 7 is uniformly bounded. Hence, 



1^43 1 < c(^-?7)-"^ 



s J^' 



X [{x - s^-'^iri - x^-^ + [x - sY-\ri - a;)"-^] ^xdS,' 



< c(i-ry)-"+-(7y~Cr 



- x)i-\x - - s)^-^dxd^' 



+c{t - 7^)-"+^T^ - C)^-^ J J {C-x)^^-^ix-e)^-\C -s)"'^dxde 

which impUes that (|6.7I ) holds for i = 3. 

Step 4 Suppose i = 4. We are going to use the following decomposition 

-{x,y)dxdy 



dxdy 



We need to compute the second derivative: 
d^G 



Hence, we can write 

d^G 



where 



ip{z) = {{a-lf -{a-l){a-2)z)(l){z) 
-(f)" (z) z^{l + z)- (j)'{z)z(l + 2z). 

We are going to use the decomposition 

ll){z) = ?Ai(z) + V'2(z), 

where 

Mz) = -{a-l)ia-2)z^{z) 

Mz) = {a-lf^{z)-(l)"{z)z\l + z)~ct>'{z)z{l + 2z). 
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This leads to 

[ A4d^di] = Bi+B2, 

•'s<i<r]<t 

where 

D 



y~x 

i = 1,2, and 

D = {{C, x,^,7],y,7]') : s <x <^ <ij <y <T]' <t}. 
Step 5 Estimation of Bi. Denote 

Di = {{^', x,tv,v') ■ s < C < X < ^ < V < V' < t}- 
Using with /3 = 1 - a - J with 5 < e/3, we obtain 

< c [ -7^r"+'-\^-a'"+''\x-sr-'{y-xr-^ 

J D 

fx — 5 \ 

X(/) d^' dxd£,dr]dydr]' 

\y ~ 

< c / (yy' - 77)-"+"-i(e - ey+'-^x - _ x)-^-^dCdxd^dT]dydT]' 
Jd 

= C f (77' - - - s)2"-2+^ 

X [{ri - x)-^-^ ~ (r]' - x)-^-^] d^dxd^dridr]' 

< c f (r]' -7]y+'-\^^^')-"+'-\x-sf''-^+\r]^x)-^-\r,-r,TdCdxd^dr^dr]' 

JDi 
Jdi 

(x - ef"'\C' - sy^^\v - - x)-'^^"^d^'dxd(dr]dr]' 

< c f (t]' - r^y-^i^ - x)-^+\x - O'^'^^^'i^' - - O'^^^d^dxdCdrjdij' 

JDi 

< 00. 

Step 6 Estimation of B2. Let us compute the function "02 (^): 
Jo 

-2q^{l + qz)-^z^{l + z) + q{l + qz)-^z(l + 2z)] dq 

\\ - + qz)-^ [{a - If {I + qzf 

+qz{l + (2 - q)z)]dq 
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This implies that the function ip2{z) is uniformly bounded. As a consequence, 
we deduce the following estimates 

\B2\ < c [ iv' -r])-''+'-\^-er''+'-\x-sr-^{y-0"-^dC'dxd^dvdydrj' 
Jd 

< c [ {r^' -y)i-\y-ri)-<-+^^-^')--+--\x-sr-^{y-^r-^d^'dxd^dr^dydr^' 

Jd 

Jd 

x(x - s)"-^irj - (,)-'^+^d^'dxd^dr]dydr]' 

< c /(r7'-y)t-i(y-r7)-i+t(e-x)t-i(x-0"'+* 



x(C' - s)-'^+i{Tj - ^)-^+id^'dxd^dr]dydr]' 

< 00. 
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